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I. INTRODUCTION 

Neutrino-nucleus cross sections have been intensively 
studied for decades. However, the precision and under- 
standing of the cross sections around 1 GcV arc still 
not satisfactory. The next generation of neutrino oscilla- 
tion experiments will search for sub-leading flavor oscilla- 
tion and charge-parity symmetry violation; the precision 
needed for these searches drives the need for more accu- 
rate independent measurements of neutrino cross sections 
[3, IH . Although several interaction channels contribute 
to the total neutrino-nucleus cross section in the neutrino 
energy range of a few GeV, an understanding of neutral 
current neutral pion production (NC7r°) is especially im- 
portant. NC7T events form the largest ^-induced back- 
ground to neutrino experiments measuring — > v e os- 
cillations in the neutrino energy range of a few GeV or 
less, such as the T2K experiment NC7r° events can 
mimic v e signal events when, for example, one of the two 
photons associated with 7r°->77 is not detected. 

NC7r°production has been measured by several past ex- 
periments 0-0] • However, their results have large uncer- 
tainty due to low statistics and are not useful expressions 
for predicting electron backgrounds in — > v e oscilla- 
tion searches, since they are typically given as ratios to 
the charged current (CC) single pion production cross 
section, which is also poorly known. T2K [l[ uses a neu- 
trino beam whose mean neutrino energy is approximately 
0.8 GeV. The experiment requires less than a 10% uncer- 
tainty on NC7T production to maintain high sensitivity 
for the Up — > v e oscillation search. Recently, two exper- 
iments published NC tt° production results. The K2K 
collaboration reported NC 7r° production in water with 
a 1.3 GcV mean neutrino energy beamQ, finding their 
measurement consistent with the Monte Carlo QV1C) pre- 
diction based on the Rein and Sehgal model Q- Mini- 
BooNE reported the yield and spectral shape of tt° s as 
a function of ir° momentum and the observation of NC 
coherent 7r° production in mineral oil (CHg) in neutrino 
beam of mean neutrino energy 0.7 GeV [lOj . The to- 
tal NC7T cross section below 1 GeV has still not been 
precisely measured yet. 

Pions are produced mainly through two distinct mech- 
anisms by neutrinos with energies around 1 GcV. In the 
dominant mode, resonant pion production, the neutrino 
interacts with a nucleon in the nucleus and excites it 
to a baryonic resonance, such as A (1232), which sub- 
sequently decays to a pion and a nucleon. The other 
mode, coherent pion production occurs when the neu- 
trino interacts with the target nucleus so that no nuclear 
breakup occurs. Resonance production and decay in a 
nuclear target differs from the case of the free nucleon 
target. This is due to nuclear effects such as Fermi mo- 
tion, Pauli blocking, and the nuclear potential. In addi- 
tion, produced mesons and baryons interact with nuclear 
matter until they escape from the target nucleus. Due 
to this final state interaction, the number, momenta, di- 
rections and charge states of produced particles can be 



changed in nuclear matter. Although there are several 
theoretical models of these processes, their uncertainties 
are still large. To understand the production mechanism 
and the nuclear effects, measurements of emitted ir° kine- 
matics are very important. 

Recent measurements of coherent pion production have 
drawn much attention. For CC coherent pion production, 
the K2K experiment placed a limit on the ratio of the CC 
coherent pi on p roduction to the total CC cross sections 
at 1.3 GeVhJ. This result was confirmed by the Sci- 
BooNE experiment [IH, although recent data from the 
SciBooNE collaboration suggest CC coherent pion pro- 
duction at a level below the pub lished limit in both neu- 
trinos and antineutrinos [l3l. Il4j. Moreover, evidence for 
NC coherent pion production with neutrino energy less 
than 2 GeV has been reported by the MiniBooNE Col- 
laboration [l(|- Hence, it is interesting to search for NC 
coherent ir° production in the SciBooNE data. 

In this paper, we present measurements of the NC7r° 
interaction in polystyrene (CsHg) using the same neu- 
trino beam as MiniBooNE (with mean neutrino energy 
0.7 GeV). We measure the ratio of the total inclusive 
NC7T cross section to the total CC cross section and kine- 
matic distributions of the emitted 7r°s. We also extract 
the fraction of coherent NC tt° events from the inclusive 
NC 7T° data sample. In these analyses, we define NC7r° 
events to be NC neutrino interactions with at least one 
7T° emitted in the final state from the target nucleus. 

II. EXPERIMENT DESCRIPTION 

A. Neutrino Beam 

The SciBooNE experiment detected neutrinos pro- 
duced by the Fermilab Booster Neutrino Beam (BNB). 
The same BNB beam is also serving the MiniBooNE 
experiment. The BNB uses protons accelerated to 8 
GeV kinetic energy by the Fermilab Booster synchrotron. 
Beam properties are monitored on a spill-by-spill basis, 
and at various locations along the BNB line. Transverse 
and directional alignment of the beam, beam width and 
angular divergence, beam intensity and losses along the 
BNB, are measured and used to monitor data quality. 
Protons strike a 71.1 cm long beryllium target, produc- 
ing a secondary beam of hadrons, mainly pions with a 
small fraction of kaons. A cylindrical horn electromag- 
net made of aluminum surrounds the beryllium target to 
sign-select and focus the secondary beam. For the data 
set used in this measurement, the horn polarity was set 
to neutrino mode, focusing (defocusing) particles with 
positive (negative) electric charge. The neutrino beam is 
mostly produced in the 50 m long decay reg ion, mainly 
from ir + — > ^ + v^ in-flight decays. See [la ] for further 
details. 
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B. Neutrino Detector 

The SciBooNE detector was located 100 m downstream 
from the beryllium target on the axis of the beam. The 
detector comprised three sub-detectors: a fully active and 
finely segmented scintillator tracker (SciBar), an electro- 
magnetic calorimeter (EC), and a muon range detector 
(MRD). SciBar served as the primary neutrino target for 
this analysis. 

Fig. [T] shows an event display of a typical NC7r° pro- 
duction event candidate. SciBooNE uses a right-handed 
Cartesian coordinate system in which the z axis is the 
beam direction and the y axis is the vertical upward di- 
rection. The origin is located on the most upstream sur- 
face of SciBar in the z dimension, and at the center of the 
SciBar scintillator plane in the x and y dimensions. Since 
each sub-detector is read out both vertically and horizon- 
tally, two views are defined: top (x vs. z projection) and 
side (y vs. z projection). 

The SciBar detector [l6[ was positioned upstream of 
the other sub-detectors. It consists of 14,336 extruded 
plastic scintillator strips which serve as the target for the 
neutrino beam as well as the active detection medium. 
Each strip has a dimension of 1.3 x 2.5 x 300 cm 3 . The 
scintillators are arranged vertically and horizontally to 
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FIG. 1. Event display of a typical NC7r° event candidate in 
SciBooNE data. The neutrino beam runs from left to right 
in this figure, encountering SciBar, the EC and MRD, in that 
order. The circles on SciBar indicate ADC hits for which the 
area of the circle is proportional to the energy deposition in 
that channel. This event display shows the electromagnetic 
shower tracks from the pair conversions of the two it decay 
photons. 



construct a 3 x 3 x 1.7 m 3 volume with a total mass 
of 15 tons. SciBar has about four radiation lengths of 
material along the beam direction. Each strip was read 
out by a wavelength shifting (WLS) fiber attached to a 
64-channel multi-anode PMT. Charge information was 
recorded for each channel, while timing information was 
recorded in groups of 32 channels by taking the logical 
OR with multi-hit TDC modules [lj- 

The gains of all PMT channels, attenuation of WLS- 
fibers, and light yield of each scintillator were continu- 
ously monitored in situ using LED and cosmic ray data 
taken between beam spills, with precision better than 
1%. The timing resolution for minimum- ionizing parti- 
cles was evaluated with cosmic ray data to be 1.6 ns. 
The average light yield for minimum-ionizing particles is 
approximately 20 photoelectrons per 1.3 cm path length, 
and the typical pedestal width is below 0.3 photoelec- 
tron. The hit finding efficiency evaluated with cosmic 
ray data is more than 99.8%. The minimum length of a 
reconstructable track is approximately 8 cm (three layers 
hit in each view). The track finding efficiency for single 
tracks of 10 cm or longer is more than 99%. 

The EC is located just downstream of SciBar, and is 
designed to measure the electron neutrino contamination 
in the beam and tag photons from ir° decay. The EC is 
a "spaghetti" type calorimeter comprised of 1 mm diam- 
eter scintillating fibers embedded in lead foil [l~8j]. The 
calorimeter is made of 64 modules of dimensions 262 x 
8x4 cm 3 . The fibers are bundled in two independent 
groups of 4 x 4 cm 2 transverse cross section, read at 
both ends by Hamamatsu PMTs. The EC comprises one 
vertical and one horizontal plane (32 modules each), cov- 
ering an active area of 2.65 x 2.65 m 2 . The EC has 
a thickness of ff radiation lengths along the beam di- 
rection. The charge information from each PMT was 
recorded. A minimum ionizing particle with a minimal 
path length deposits approximately 91 MeV in the EC. 
The energy resolution for electrons was measured to be 
14%/y/E (GcV) using a test beam [ij]. The detection 
efficiency for cosmic muons is 96%; the inefficiency stems 
from gaps between the modules. 

The MRD was installed downstream of the EC and 
is designed to measure the momentum of muons pro- 
duced by CC neutrino interactions. It had 12 iron plates 
with thickness 5 cm sandwiched between planes of 6 mm 
thick scintillation counters; there were 13 alternating hor- 
izontal and vertical planes read out via 362 individual 
2 inch PMTs. Each iron plate measured 274 x 305 cm 2 . 
The MRD measured the momentum of muons up to 
1.2 GeV/c using the observed muon range. Charge and 
timing information from each PMT were recorded. Hit 
finding efficiency was continuously monitored using cos- 
mic ray data taken between beam spills; the average hit 
finding efficiency is 99%. 
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C. Data Summary 

The SciBooNE experiment took data from June 2007 
until August 2008. After applying data quality cuts to 
all beam events (l2|. 2.52 x 10 20 protons on target are 
usable for physics analysis, with 0.99 x 10 20 protons on 
target collected in neutrino mode. The analysis presented 
herein uses only neutrino mode data. 

III. EXPERIMENT SIMULATIONS 

A. Neutrino Flux Prediction 

Predictions for the BNB neutrino flux illuminating the 
SciBooNE detector are obtained via a GEANT4 simu- 
lation of the beamline. The simulation accounts for all 
relevant beamline geometry and materials, the measured 
BNB beam optics properties, and the horn magnetic field. 
Hadronic interactions are carefully modeled. Cross sec- 
tions for elastic, quasi-elastic and other inelastic interac- 
tions of charged pions and nucleons with beryllium and 
aluminum are simulated according to a custom model 
validated with external data, and covering the most rel- 
evant momentum range (down to 0.5 GeV/c for pions, 
2 GeV/c for nucleons). Furthermore, the multiplicity 
and kinematics of all relevant particle types produced 
in the inelastic interactions of primary (8.4-8.9 GeV/c) 
protons with beryllium are also described by a custom 
model based on external data. For 7r + production, a pa- 
rameterization based on HARP and BNL E910 |20l| 
data is used. Other hadronic and all electromagnetic 
processes of importance to neutrino production are de- 
scribed instead by standard GEANT4 models. The mod- 
eling of neutrino-producing weak decays incorporates ac- 
curate knowledge of meson decay branching fractions and 
form factors, and includes muon polarization effects. For 
a detailed description of the BNB simulation code, see 
[l5j |. According to the simulation, the neutrino flux at 
the SciBooNE detector is dominated by muon neutrinos 
(93%), while the neutrino energy spectrum peaks at ~ 
0.6 GeV, has a mean neutrino energy of ~ 0.7 GeV, and 
extends up to 2-3 GeV [lsj . 

B. Neutrino Interaction Simulation 

In the SciBooNE experiment, neutrino interactions 
with nuclear targ ets are simulated by the NEUT pro- 
gram library [2l[ 1221] that is used in the Kamiokandc, 
Super-Kamiokande, K2K, and T2K experiments. 

The nuclear targets handled in NEUT are protons, car- 
bon, oxygen, and iron. The energy of neutrinos handled 
by the simulation ranges from 100 MeV to 100 TeV. 
The types of neutrino interactions simulated in both 
NC and CC are : elastic and quasi-elastic scattering 
(vN — > £N'), single meson production (vN — ► £N'm), 



single gamma production {vN — > £N'j), coherent ir pro- 
duction (i/ 12 C(or 56 Fe) -> £tx 12 C(or 56 Fe)), and deep in- 
clastic scattering (vN — > iN'hadrons), where N and 
N' are the nucleons (proton or neutron), £ is the lep- 
ton (electron, muon or neutrino), and m is the meson. 
In nuclei, interactions of the mesons and hadrons with 
the nuclear medium are simulated following the neutrino 
interactions. 



1. Single meson production via baryon resonances 

The main signal in this analysis is NC single ir° produc- 
tion via baryon resonances. The resonant single meson 
production is simulated based on the model of Rein and 
Sehgal H- The cross section of the NC single ir° pro- 
duction per nucleon on a polystyrene target (CgHs) in 
NEUT is shown in Fig. [2] The per nucleon cross section 
of a CH molecule is calculated by summing the contribu- 
tions from the six protons and six neutrons bound in the 
carbon nucleus as well as the free proton, and dividing 
that by 13. Following production, the intra-nuclear in- 
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FIG. 2. Cross section versus neutrino energy of NC single 
7T° production per nucleon on a polystyrene target (CsHg) 
estimated in NEUT; the curve is based on the Rein and Sehgal 
model. 

teractions of the meson and nucleons are simulated using 
a cascade model in which the particles are traced until 
escaping from the nucleus. According to this model, ap- 
proximately 40% of 7r°s interact in the target nucleus, 
averaged over our neutrino flux. For scattering off nucle- 
ons in the nucleus by a neutrino, the relativistic Fermi 
gas model of Smith and Moniz [23[ is implemented. The 
nucleons are treated as quasi-free particles and the Fermi 
motion of nucleons along with the Pauli exclusion princi- 
ple is taken into account. The Fermi surface momentum 
is set to 217 MeV/c and the nuclear potential is set to 
27 MeV for carbon. The vector and axial- vector form fac- 
tors are formalized to be dipole with Mj." = 0.84 GeV/c 2 
and Mjj ir = 1.21 GeV/c 2 . The same Fermi momen- 
tum distribution, nuclear potential and Q 2 -dependencc 
of form factors are used in all other neutrino-nucleus in- 
teractions except for coherent ir production. 

The Rein and Sehgal model assumes an intermedi- 
ate baryon resonance, N* , in the reaction of vN — > 
£N*,N* —> N'm. All intermediate baryon resonances 
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with mass less than 2 GeV/c 2 are included. Baryon reso- 
nances with mass greater than 2 GeV/c 2 are simulated as 
deep inelastic scattering. Pion-less A decays — which pro- 
duce nopion in the final state and account for 20% of A 
events [241 ] — are also simulated. To determine the angu- 
lar distribution of final state pions, Rein's method [25| is 
used for the P33(1232) resonance. For other resonances, 
the directional distribution of the generated pion is cho- 
sen to be isotropic in the resonance rest frame. 

The inelastic scattering, charge exchange and absorp- 
tion of pions in nuclei are simulated. The interaction 
cross sections of pions in the nuclei are based on the 
model by Salcedo et al. [26j . For inelastic scattering and 
charge exchange interactions, the direction and momen- 
tum of pions are affected. In the scattering amplitude, 
Pauli blocking is also taken into account. 

2. Coherent n production 



TABLE I. The expected number and fraction of events in each 
neutrino interaction estimated by NEUT at the SciBooNE de- 
tector location with the neutrino beam exposure of 0.99 x 10 20 
protons on target. The 10.6 ton fiducial volume of the SciBar 
detector is assumed. CC and NC interactions are abbreviated 
as CC and NC, respectively. 



Interaction Type # Events Fraction(%) 



CC quasi-elastic 


53,363 


41.4 


CC single it via resonances 


29,688 


23.1 


CC coherent n 


1,771 


1.4 


CC single meson except n 


839 


0.7 


CC DIS 


6,074 


4.7 


NC elastic 


22,521 


17.5 


NC single 7r° via resonances 


6,939 


5.4 


NC coherent n° 


1,109 


0.9 


NC single meson except iy° 


4,716 


3.7 


NC DIS 


1,768 


1.4 



The 7T° signal events contain a contribution from NC 
coherent ir° production. Due to the small momentum 
transfer to the target nucleus, the outgoing neutrino and 
the pion tend to go in the forward direction. The formal- 
ism developed by Rein and Sehgal [13, 28] is used to sim- 
ulate the interactions. The axial vector mass, M'^' herent , 
is set to 1.0 GeV/c 2 , and the nuclear radius parameter R 
is set to 1.0 fm. For the total and inelastic pion-nucleon 
cross sections used in the formalism, the fitted results 
given in Rein and Sehgal's paper are employed. The NC 
coherent tt° production cross section on a polystyrene tar- 
get is shown in Fig. with the NC single ir° production 
via baryon resonances and the total CC cross sections. 
The Rein and Sehgal model predicts the NC coherent 
7T° production rate to be approximately 1% of the total 
neutrino CC rate in SciBooNE. 
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3. Quasi-elastic scattering and deep inelastic scattering 

The dominant interaction in the SciBooNE neutrino 
energy is CC quasi-elastic scattering, which is imple- 
mented using the Smith and Moniz model |23j]. My E 
and Ml E are set to be 0.84 GcV/c 2 and 1.21 GeV/c 2 , 
respectively. 

The deep inelastic scattering (DIS) cross section is 
calculated using the GRV98 parton distribution func- 
tions [29| . As well as quasi-elastic scattering and deep 
inelastic scattering, other neutrino interactions in NEUT 
are described in [12( in detail. 

With the SciBooNE neutrino beam exposure of 0.99 x 
10 20 protons on target, the expected number of events 
in the SciBooNE detector in each neutrino interaction 
is listed in Tab. |U For the purpose of systematic stud- 
ies of neutrino interaction simulations, we also use the 
NUANCE event generator [30| that is used in the Mini- 
BooNE experiment. The types and models of neutrino 
interactions in NUANCE are similar to those of NEUT 
but with different treatment of re-interactions of mesons 
and hadrons with the nuclear medium. 



C. Neutrino detector simulation 



FIG. 3. Cross sections versus neutrino energy of the to- 
tal CC interaction (solid line), the NC single tv° production 
via baryon resonances (dashed line) and the NC coherent 7T° 
production (dotted line) calculated in NEUT per nucleon on 
polystyrene target. 



The GEANT4 framework is used for the detector simu- 
lation. The Bertini cascade model within GEANT4 [U 
is used to simulate the interactions of hadronic parti- 
cles with detector materials. The detector simulation 
includes a detailed geometric model of the detector, in- 
cluding the detector frame and experimental hall and 
soil, based on survey measurements taken during detec- 
tor construction. A description of the detector simulation 
is given in (i~2l ]. 

In addition to neutrino interactions inside the detector, 
we also simulate interactions in the surrounding material 
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(the walls of the detector hall and soil). The density of 
material is assumed to be 2.15 g/cm 3 for the calculation 
of the interaction rate, and concrete of that density is 
used as the material for propagation of product particles. 
We generate events in a volume of ±5 m in x, y, and z 
direction in the SciBooNE coordinates. 



IV. DATA ANALYSIS 

The present analysis has two main goals. The first is 
to measure the ratio of NC7r° production cross section to 
the total CC cross section. We measure the ratio of cross 
sections in order to minimize systematic uncertainty due 
to the neutrino flux prediction. The second goal is to 
measure the ir° momentum and angular spectra. In ad- 
dition to the two main goals, we also extract the coherent 
7T° fraction in the context of the Rein and Sehgal model. 

We reconstruct gamma rays converting in SciBar 
and select events with two reconstructed gamma rays 
and no muons, which is the characteristic topology of 
NC7r°events. We do not include NC7r°events in which 
one or both gamma rays convert in the EC. 



A. Signal Definition 



We define an NC7r° interaction as an NC neutrino in- 
teraction in which at least one ir° is emitted in the final 
state from the target nucleus, v^C — > v^X where X 
represents the nuclear remnant and any combination of 
nucleons and mesons. According to our MC simulation, 
96% of NC7T events without any selection cuts have a sin- 
gle 7T° (85% from a single ir° without any other mesons 
and 11% from a single ir° with charged mesons) and 4% 
have two 7r°s. Any ir° emitted from the initial target nu- 
cleus constitutes a signal event whether it is created from 
the neutrino vertex or final state interactions. Events 
with a 7T° produced in the neutrino interaction but ab- 
sorbed in the target nucleus are not included in the signal 
sample, nor are events in which 7r°s are produced by sec- 
ondary particles interacting with the detector scintillator 
outside the target nucleus. 



B. Gamma Ray Reconstruction 

1. Gamma Conversion Probability 

Since the length of SciBar in the beam direction cor- 
responds to four radiation lengths, a significant fraction 
of gamma rays escape from SciBar without conversion. 
In 30% of events with a 7r° emitted within SciBar's fidu- 
cial volume, both gamma rays convert in SciBar; in 38%, 
only one gamma ray converts in SciBar; in 32%, neither 
gamma ray converts in SciBar. Since we aim to recon- 
struct two gamma rays to identify NC7r° events, the max- 
imum detection efficiency attainable is 30%. 



2. Track Reconstruction 

The first step of the event reconstruction is to search 
for two-dimensional tracks in each view of SciBar using a 
cellular automaton algorithm 32j. For tracking, the hit 
threshold is set to 2.5 photo-electrons, corresponding to 
roughly 0.25 MeV. Three dimensional tracks are recon- 
structed by matching the timing and z-edges of the two 
dimensional projections. In order to match track projec- 
tions in a three dimensional track, the timing difference 
between two two-dimensional projections is required to 
be less than 50 ns, and the z-edge difference must be less 
than 6.6 cm for upstream and downstream edges. This 
method is used for all charged particles and is the first 
step of gamma ray reconstruction. 



3. Particle Identification Parameter 

The SciBar detector has the capability to distinguish 
protons from other particles using dE/dx since recoil 
protons at SciBooNE energies interact well above min- 
imum ionizing energy deposit. We define a muon con- 
fidence level (MuCL) using the observed energy deposit 
per layer for all reconstructed tracks [lH. The MuCL of 
the proton tracks tends to be close to while the MuCL 
of other tracks tends to be close to 1. Proton-like tracks 
arc- defined to have MuCL less than 0.03. 



4- Extended Track 

Single reconstructed tracks are extended in two ways to 
improve the energy reconstruction of gamma rays within 
SciBar. The first step is merging two or more tracks if 
they are nearly co-linear, because electromagnetic show- 
ers can form separate hit clusters in SciBar resulting in 
two or more tracks. The second step is collecting lone hits 
around merged tracks. Electromagnetic showers some- 
times deposit energy around the main track and these 
hits are missed by the track reconstruction algorithm. 
Hits not assigned to any track within 20 cm from 2D 
projections of the merged track (i.e., after the first step) 
for each view are added to the extended track. The meth- 
ods described above are applied only to non-proton-likc 
tracks. For the energy reconstruction, we use charge in- 
formation of hits associated with original tracks as well as 
newly assigned hits. For reconstructing the directions of 
gamma rays, we fit positions of hits in all original tracks 
in the extended track with a straight line and do not use 
hits newly collected in the second step. 



5. SciBar-EC matching 

Gamma rays can escape SciBar to deposit energy in 
the EC. After event reconstruction in SciBar, we search 
for EC clusters aligned with tracks from SciBar. One EC 
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cluster is denned as a collection of neighboring EC hits. 
For an EC hit, the pulse heights of both side PMTs are 
required to be above threshold, which is set to 3 times 
the width of each pedestal — about 7 MeV. The energy of 
an EC hit is the geometric average of the two PMTs. The 
center of an EC cluster is defined as the energy- weighted 
average of hits in the cluster. To match an EC cluster to 
an SciBar track, the EC cluster is required to be within 
10 cm of the extrapolated two-dimensional projections 
of the SciBar track in each EC plane. The energy of 
matched EC clusters is added to the corresponding ex- 
tended tracks. 



6. Gamma Ray Reconstruction Performance 



muon track on the most upstream layer of SciBar to elim- 
inate neutrino-induced incoming particles from the up- 
stream wall or soil. The neutrino interaction vertex for 
CC events is reconstructed as the upstream edge of the 
muon track. We select events whose vertices are in the 
SciBar fiducial volume, defined to be ±130 cm in both 
the x and y dimensions, and 2.62 cm< z < 157.2 cm, a 
fiducial mass of 10.6 tons. The time of the muon track 
is required to be within a 2 fis window around the beam 
pulse. Finally, we require the muon track to stop in the 
MRD. The MRD-stopped event sample [l2[ serves as the 
normalization sample in the cross-section ratio measure- 
ment. Unless otherwise indicated, the MC distributions 
in this paper are normalized using the MRD-stopped 
data sample. 



We study the performance of the gamma ray recon- 
struction algorithms using the MC simulation. The 
angular resolution of gamma rays passing all selection 
cuts (see section HV D[) is estimated to be approximately 
6°. For the energy reconstruction, energy of matched 
EC clusters is added to corresponding extended tracks. 
About 7% of selected extended tracks are made by two 
or more tracks. The average gamma ray energy deposit 
of SciBar is estimated to be 116 MeV and the energy 
resolution is estimated to be 6%. About 17% of selected 
extended tracks have matched EC clusters. The average 
gamma ray energy deposit in such matched EC clusters 
is estimated to be 72 MeV and the energy resolution is 
estimated to be 32%. 

Not all gamma ray energy is deposited in, nor recorded 
by, the detector. Such lost energy is called leakage. The 
actual leakage, defined as 

gamma energy in extended track 

-L a ct = 1 — 3 (1) 

true gamma energy 

is estimated to be 24%; 11% comes from energy loss in 
passive regions and gamma rays escaping from the de- 
tectors and 13% comes from energy deposit in active re- 
gions but not assigned to the extended track. The recon- 
structed energy can include energy deposited by other 
particles. On average, 15% of the total energy in an ex- 
tended track comes from other particles. Due to this 
contamination, the effective leakage, defined as: 

reconstructed energy of extended track , , 

L eff = 1 7 ~ (2) 

true gamma energy 

is 15%, which is smaller than L act averaged over all NC7r° 
events. 



C. CC Event Selection 

To identify CC events, we search for events in which 
at least one reconstructed track in SciBar, when pro- 
jected out of SciBar, is matched with a track or hits in 
the MRD. We reject events with hits associated with the 



D. NCtt" Event Selection 

The clearest feature of NC ir° production is two gamma 
rays, coming from the decay of the 7r°, converted into two 
e + e~ pairs. Background events stem from sources both 
internal and external to SciBar. Internal backgrounds 
are neutrino interactions other than NC7r° (mainly CC) 
within SciBar; external backgrounds come from neutrino 
interactions in the material outside of the detector vol- 
ume (support structure, walls and soil — so called dirt 
events) as well as from cosmic rays. In dirt events, neu- 
trinos interact with materials such as the walls of the 
experimental hall or soil and produce secondary parti- 
cles which deposit energy within SciBar's fiducial volume 
and cause recorded hits. The contribution of accidental 
cosmic rays in any event sample is small and accurately 
estimated by data taken with off-beam timing; the frac- 
tion of accidental cosmic ray events is 1.8% after all event 
selection cuts. Data distributions shown hereafter have 
been cosmic-ray subtracted from them. The event selec- 
tion cuts for NC7T production are developed to select 
events with two gamma rays while rejecting these back- 
grounds. 

1. Preselection Cuts 

We select events with at least two three-dimensional 
tracks. The timing of tracks are required to be within 
50 ns of each other and the closest endpoints of any 
two tracks are required to be in the fiducial volume de- 
fined in Sec. IIVD1 The closest endpoints are chosen as 
the edge combination for two given tracks that gives the 
minimum distance among the four possible combinations. 
The times of the two tracks are required to be within the 
2 /zs beam window. In addition, we reject events with 
hits in the first layer of SciBar within 100 ns from the 
times of tracks, to remove dirt background events. 
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FIG. 4. Atmax distribution (|IVD2[) after rejection of side 
escaping tracks. The contributions from events with decay 
electrons, events without decay electrons and the dirt events 
are shown separately for the MC simulation. 

2. Muon Track Rejection Cuts 

Events with rmions are predominantly CC events. To 
avoid muons which escape SciBar but do not penetrate 
the MRD, we reject events with tracks escaping from the 
side of SciBar; both edges of tracks must be in ±130 cm 
in both the x and y dimensions. 

To reject muons stopping in SciBar, we tag the elec- 
trons from muon decay. No charge information is 
recorded for any scintillator strip after its first hit in an 
event, but the times of hits above threshold continue to 
be recorded. Thus most decay electrons are not recon- 
structed as tracks but can be identified as delayed time 
hits near the end of a muon track. We examine the max- 
imum timing difference (Ai max ) between the track times 
and late hits times at the ends of tracks. Events with 
decay electrons yield high values of At max because of the 
long muon lifetime (t m = 2.2 /i sec). Figure 2] shows 
the At max distribution. We select events with Ai max less 
than 100 ns. 



FIG. 5. The minimum distance between two tracks after 
muon rejection cuts. The contribution from NCw produc- 
tion, internal backgrounds with a it in the final state, in- 
ternal background without a 7r° in the final state and "dirt" 
background events are shown separately for the MC simula- 
tion. 

4- Electron Catcher Cut 

Matched EC Clusters are used to reject muons pen- 
etrating the EC. Two quantities are used: the energy 
deposit in matched EC clusters in the upstream layer, 
Ei (E 2 is energy in the downstream layer), and the ra- 
tio of energy deposits in the downstream EC cluster over 
the upstream EC cluster, Rec = E2/E1. Figure O [7] and 
[8] shows the distributions of the number of tracks with 
matched EC clusters, E\ and Rec, respectively. Events 
satisfying one of the three following condition are se- 
lected, (i) No matched EC clusters or (ii) E x > 150 MeV 
or (hi) Rec < 0.2. For events without matched EC clus- 
ters, both Ei and Rec are left undefined (not included 
in Fig. [7] and [5]). For events with only upstream matched 
EC clusters, the minimum E\ is chosen and Rec is sc t 
to 0. For events with both upstream and downstream 
matched EC clusters, the maximum Rec is chosen and 
the corresponding E\ is chosen. 



3. Track Disconnection Cut 

CC events often have two or more tracks with a com- 
mon vertex while two gamma rays from 7r°s are typically 
isolated from each other. Hence, the distance between 
two tracks is used to separate two gamma rays from CC 
events. If there are two particles with a common vertex, 
the minimum distance is close to zero. Figure [5] shows 
the distribution of the minimum distance. Events with a 
minimum distance greater than 6 cm are selected. 



5. Number of Photon Tracks 

We use the extended track information instead of the 
track information hereafter. As described in Scc lIV BT1 
we only use non-proton-like tracks to obtain extended 
tracks. Fig. |H] shows the MuCL distribution after the 
EC cut. The dashed line (MuCL = 0.03) separates into 
particles to proton-like or non-proton-like tracks. The 
gamma ray efficiency for the non-proton-like sample is 
87% and the contamination of gamma rays in the non- 
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FIG. 6. Option 1 of the EC cut: the number of track with 
matched EC clusters after the track disconnection cut. Events 
without SciBar-EC matched tracks pass the EC cut. Events 
with SciBar tracks matching EC clusters can pass the EC cut 
if they satisfy the requirement shown in Fig. [7] or [8] 




400 
Ei (MeV) 



FIG. 7. Option 2 of the EC cut: the energy deposited in the 
upstream layer of the EC (Ei) after the track disconnection 
cut. Events with E\ > 150 MeV pass the EC cut. Events 
with Ei < 150 MeV can pass the EC cut if they satisfy the 
requirement shown in Fig. H 



proton-like sample is 81%. 

Figure [TUJ shows the distribution of the number of ex- 
tended tracks after the EC cut. To reconstruct 7r°s, 
events with more than one extended track are selected. 
As shown in Fig. [TUJ 58% of the NC7r° events have only 
one extended track and are rejected by this cut. How- 
ever, 39% are events with one reconstructed gamma ray, 



FIG. 8. Option 3 of the EC cut: the ratio of deposited energy 
in the downstream over the upstream layer (R) after the track 
disconnection cut. Only events with Ei < 150 MeV are 
shown. 




FIG. 9. The muon confidence level (MuCL) of tracks after the 
EC cut. The contributions from true gamma ray, muon, and 
proton tracks are shown separately for the MC simulation. 



due to mis-reconstruction or gamma rays not converting 
in SciBar. Such events can not be used for 7r° recon- 
struction; 12% are events where the two gamma rays are 
reconstructed as two tracks but one of them is identified 
as a proton-like track; 7% are events in which the two 
gamma rays are reconstructed as two tracks but then 
merged, resulting in one extended track. 

This cut is also effective at rejecting dirt backgrounds 
since many dirt background events have only one ex- 
tended track, as shown in Fig. [TU] 
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FIG. 10. The distribution of the number of extended tracks 
after the EC cut. 



6. Reconstructed ir° Vertex Position Cut 

The reconstructed vertex position of a ir° is calcu- 
lated as the intersection of two two-dimensional extended 
tracks. First, we calculate the intersection point for all 
combinations of extended tracks in each view — two z po- 
sitions (ztop and Zside) for each combination. We choose 
the combination giving the minimum | z top — z S id c | as the 
7r candidate. The reconstructed z-vertices are obtained 
by taking the error weighted average of Zt op and z s ide^ 



Ztop 




Sz top 2 




1 


1 







(3) 



where 5z top ( S ide) is the error on z top ( sido ) returned by the 
track reconstruction algorithm. Figure [TTJ shows the re- 
constructed z-vertices of 7r°s. The vertex resolution is 
approximately 12 cm for all three dimensions. Most 
events with a ir° produced in SciBar yield a vertex within 
SciBar — but many dirt events yield a vertex position up- 
stream of SciBar — so we select events with reconstructed 
7T° z-vertex greater than cm. 



7. Reconstructed n° Mass 

Figure[T2]shows the reconstructed mass of the ir° calcu- 
lated as y^E^E^l -cos6"' cc ), where E^f and E'f are 

the reconstructed energies of the extended tracks (E™ > 
E™ 2 C ) and 9 lec is the reconstructed angle between the ex- 
tended tracks. The MC simulation describes well the tail 
of the distribution, which is background-dominated. We 
select events with 50 MeV/c 2 < M™ C < 200 MeV/c 2 . 
The peak value is smaller than the actual 7r° mass (135 
McV) due to energy leakage of 7s. 
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FIG. 11. The reconstructed z-vertices of 7r°s after the require- 
ment of at least two extended tracks. 
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FIG. 12. The reconstructed mass of 7r°s after the recon- 
structed vertex position cut. 



8. Event Selection Summary 

Table QTJ shows the number of events in data and the 
MC at each event selection stage. The numbers for the 
MC simulation arc normalized to the number of MRD 
stopped events. We select 657 events after all cuts. Sub- 
tracting the estimated background of 240 events (202 in- 
ternal and 38 external) yields 417 signal events. The MC 
expectation is 368 events. The purity of NC ir° produc- 
tion after all event selection cuts is estimated to be 61% 
(40% from single n production via resonance decay, 15% 
from coherent n production and 5% from neutrino deep 
inelastic scattering). According to our MC simulation, 
96% of selected NC71" events have one ir° (91% from a 
single 7T° without any other mesons and 5% from a sin- 
gle 7T° with charged mesons) and 4% have two ir°s. The 
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efficiency for NC7r° production, defined as: 

the number of selected NC7r° events 



e NCir° 



the number of generated NC7r° events 



, (4) 



is estimated to be 5.3%. The internal background, which 
accounts for 33% of this sample, contains CC 7r° produc- 
tion including secondary ir°s (18%), NC secondary ir° 
production in detector materials (9%) and non-7r° back- 
ground (6%). According to our MC simulation, the av- 
erage energy of neutrinos producing NC7r° events in the 
SciBar fiducial volume is 1.3 GeV and the average en- 
ergy of neutrinos producing NC7r° events that pass all 
selection cuts is 1.1 GeV. The average energy of neutri- 
nos producing NC7r°s coherently is 1.1 GeV, while the 
average energy that pass the selection cuts is 1.0 GeV. 



V. RESULTS 

A. cr(NC7r°)/cr(CC) Cross Secion Ratio 

We measure the ratio of the NC71" production to the 
total CC interaction cross sections. 



1. NCtt Production 

The efficiency corrected number of NC7r° events is cal- 
culated as: 



7V(NC7r°) = 



iVobs - N BG 

£ NC7T 



(5) 



where N b s is the number of observed events, Nbg is 
the number of background events estimated by the MC 
simulation, and e NCjr o is the selection efficiency of NC 
7T° events calculated by the MC simulation. Nobs and 
Nbg, £ NC7r° are 657, 240.0 and 0.053, respectively. After 
subtracting background and correcting for the selection 
efficiency, the number of NC7r° candidates is measured 
to be [7.8 ± 0.5(stat.)] x 10 2 . For the background cal- 
culation, we use the MC expectation normalized to the 
number of MRD stopped events. The neutrino energy 
dependence of the selection efficiency for NC7r° events is 
shown in Fig. 1131 The mean neutrino energy for NC7r° 
events in the sample is estimated to be 1.1 GeV after 
event selection cuts. 



2. Total CC Interactions 

The total number of CC interactions is estimated us- 
ing the MRD stopped sample. The mean neutrino energy 
of MRD-stopped events is estimated to be 1.2 GeV. The 
number of CC candidates after correcting for the selec- 
tion efficiency is calculated as: 



N(CC) 



*S2 



ecc 



(6) 



where is the number of observed CC events, N^q 

and ecc are the number of background events and selec- 
tion efficiency in the sample, respectively, estimated with 
the MC simulation. We observed 21,702 MRD-stopped 
events (N^). The number of background events and 
the selection efficiency are estimated to be 2348 (N^q) 
and 19% (ecc), respectively. The neutrino energy depen- 
dence of the selection efficiency for CC events is shown 
in Fig. [14] After subtracting the background events and 
correcting for the efficiency, the number of CC events is 
measured to be [1.02 ± 0.01(stat.)] x 10 5 . 



3. Cross section ratio 

The ratio of the NC 7r° production to the total CC 
cross section is measured to be: 

cr(NC7r°) V(NCtt ) 
ct(CC) ~ JV(CC) 

= (7.7 ± 0.5(stat.) ± 0.5(sys.)) x 10" 2 , (7) 

at the mean neutrino energy of 1.14 GeV; systematic un- 
certainties are described in section fV A 41 The MC expec- 
tation based on the Rein and Sehgal model is 6.8 x 10~ 2 . 
The total uncertainty, ±0.7 x 10~ 2 is obtained adding 
statistical and systematic uncertainties in quadrature. 
Although the value of this measurement is larger than 
the expectation by 11%, the excess corresponds to 1.3 
standard deviations if the total uncertainty is taken into 
account. 



4- Systematic errors 

The sources of systematic error are divided into four 
categories, (i) detector response and track reconstruction, 
(ii) nuclear effects and neutrino interaction models, (iii) 
neutrino beam and (iv) dirt background. We vary these 
sources within their uncertainties and take the resulting 
change in the cross section ratio as the systematic un- 
certainty of the measurement. Table [TTT1 summarizes the 
systematic errors in the NC7r° cross section ratio. The 
total systematic error is ±0.5 x 10~ 2 on the cross section 
ratio. 

a. Detector response and track reconstruction The 
crosstalk of the MA-PMT was measured to be 3.15±0.4% 
for adjacent channels and is varied within the measure- 
ment error. The single photoelectron resolution of the 
MA-PMT is set to 50% in the simulation, to repro- 
duce the observed dE/dx distribution of cosmic muons. 
The absolute error is estimated to be ±20%. Hence, we 
vary the single photoelectron resolution by ±20%. Birk's 
constant for the SciBar scintillator was measured to be 
0.0208 ± 0.0023 cm/MeV [H and is varied within the 
measurement error. The hit threshold for track recon- 
struction is varied by ±20%. A 10% difference of the total 
pion-carbon cross section is seen for higher energy pions 
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TABLE II. Event selection summary; for the MC expectation, NC7r°signals, integral backgrounds (BG) and dirt backgrounds 
are shown separately. The number of CC background events is shown in parentheses. 



Event selection 


DAI A 




MC 




NCtt u 






NCtt 


signal Internal BG (CC) Dirt BG Efficiency 


Purity 


Pre-selection Cuts 


11,926 


1,893 


9,808 (9050) 


895 


27.3% 


15% 


Muon Track Rejection Cuts 


5,609 


1,377 


3,785 (3326) 


606 


19.8% 


24% 


Track Disconnection Cuts 


3,614 


1,314 


1,706 (1306) 


595 


18.9% 


36% 


Electron Catcher cut 


2791 


1202 


1088 (714) 


579 


17.3% 


42% 


Number of Photon Tracks 


973 


443 


389 (294) 


121 


6.5% 


46% 


7t° Reconstructed ir° Position Cut 


905 


428 


382 (288) 


65 


6.2% 


49% 


Reconstructed 7r° mass 


657 


368 


202 (140) 


38 


5.3% 


61% 



=3 

CO 
1— 

<: 



NCtt generated in FV 



NCtt selected 




4000 
Ev (MeV) 



o 
c 

o 

it 

LU 



o.oer 




4000 
Ev (MeV) 



c 
CD 




4000 
Ev (MeV) 



>> 
o 
c 

Q) 
O 
it 
LU 



0.2 



0.1" 



2000 



4000 
Ev (MeV) 



FIG. 13. Expected neutrino energy spectra and selection ef- 
ficiency as a function of neutrino energy for NCw events. 



FIG. 14. Expected neutrino energy spectra and selection ef- 
ficiency as a function of neutrino energy for all CC events. 
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TABLE III. Summary of the systematic errors in the NCw 
cross section ratio. 



Source 


error 


(xicr*) 


Detector response 


-0.39 


0.38 


v interaction 


-0.25 


0.30 


Dirt background 


-0.10 


0.10 


v beam 


-0.11 


0.22 


Total 


-0.48 


0.54 



between the GEANT4 simulation and external measure- 
ments. Hence, we vary the cross section by ±10%. The 
uncertainty of the energy scale of gamma rays is esti- 
mated to be ±3%. We vary the reconstructed energy of 
extended tracks by ±3%. For the uncertainty on recon- 
struction of the gamma direction, we study how the dif- 
ference between data and MC distributions changes when 
we change the gamma direction reconstruction algorithm. 
We compare our standard algorithm with gamma di- 
rection reconstruction obtained using extended tracks 
- resulting in poorer angular resolution. We take this 
change as the uncertainty. The largest contribution to 
the uncertainty in the cross section ratio are the crosstalk 
of the MA-PMT (-0.00.+0.27) and the hit threshold (- 
0.25,+0.05). 

b. Neutrino interaction models and nuclear effects 
The uncertainty in CC resonant pion production is es- 
timated to be approximately ±20% based on the K2K 



measurement [34|. We vary the cross section of CC res- 
onant pion production by ±20% and take that change 
as the systematic error. We also vary the NC/CC ra- 
tio by ±20% and take that change as a systematic error. 
The uncertainty in the axial vector mass for CC quasi- 
clastic and NC elastic scattering as well as CC(NC) res- 
onant pion production is estimated to be approximately 
±0.1 GeV/c 2 based on recent measurements [H, HH; re- 
sults from past experiments are systematically lower than 
recent measurements j37| . and thus we only vary Ma 
down to 1.11 GeV/c 2 , and take that change as the sys- 
tematic error. The biggest contribution to the uncer- 
tainty of the cross section ratio is the CC resonant pion 
production (-0.14.+0.16). 

We consider uncertainties in the pion interactions in- 
side the nucleus. For pions produced by neutrino interac- 
tions, uncertainties on the cross sections for pion absorp- 
tion, pion inelastic scattering and pion charge exchange 
in the nucleus are approximately 30% HH in the momen- 
tum range of pions from A decays; we vary these pion 
interaction cross sections and take the resultant change 
in the ratio as the uncertainty. The largest contribution 
to the uncertainty of the cross section ratio is the pion 
absorption (-0.17,+0.19). 

As a cross check, we measure the cross section ratio 
using the NUANCE event generator [30[ to predict event 
rates and calculate efficiencies, and obtain a measured 
ratio of 7.9 x 10~ 2 (the NUANCE expectation is 7.1 x 



10" 2 ). The result using NUANCE agrees with the NEUT 
result (7.7 x 10~ 2 ) within the systematic uncertainty, so 
we do not add the NEUT/NUANCE difference to the 
systematic uncertainty. 

c. Dirt Backgrounds As shown in Fig. I X 1 1 the dirt 
background simulation describes data at z < — 20 cm 
where the dirt background is the dominant contamina- 
tion. However, the statistical uncertainty is large, 15%. 
We scale the dirt contamination by ±15% in the final 
sample and take the change as the systematic error due 
to dirt backgrounds. 

d. Neutrino Beam The uncertainties in secondary 
particle production cross sections in proton-beryllium in- 
teractions, hadronic interactions in the target or horn, 
and the horn magnetic field model arc varied within their 
externally estimated error bands. Detailed descriptions 
of each uncertainty are found elsewhere fl5| . Systematic 
uncertainties in the neutrino flux are reduced by remov- 
ing the model dependent parameterization in the propa- 
gation of errors from the HARP data 15j . Uncertainties 
associated with the delivery of the primary proton beam 
to the beryllium target and the primary beam optics, 
which result in an overall normalization uncertainty, are 
not considered in this analysis since they cancel in the 
cross section ratio. 



B. Reconstructed ty° Kinematics 

After all event selection cuts, we studied the recon- 
structed kinematics of the 7r°s: the ir° momentum and 
cosine of the ir° angle with respect to the beam direction, 
as shown in Fig. [15] and Fig. [16] The NC7r° efficiency as 
functions of ir° momentum and angle are shown in Fig. 1171 
and 1181 respectively. The average momentum of recon- 
structed 7r°s is estimated to be 223 MeV/c while the av- 
erage momentum of true 7r°s after all event selection cuts 
is 264 McV/c according to our MC simulation. This dif- 
ference comes from energy leakage of gamma rays. The 
relation between the true and reconstructed ir° momen- 
tum is shown in Fig. 1191 The momentum resolution is 
estimated to be 23%. The relation between the true and 
reconstructed tt° direction is shown in Fig. [50] The angu- 
lar resolution of 7r°s is estimated to be 6°. In events with 
two 7r°s, we choose the ir° with the largest momentum 
when comparing the true and reconstructed kinematic 
quantities. 

In the reconstructed 7r° momentum and angular dis- 
tributions, we extract NC7r° signal events by subtract- 
ing the expected backgrounds; to estimate backgrounds, 
we use the MC expectation normalized to the number of 
MRD-stopped events. After the background subtraction, 
we convert the reconstructed ir° momentum (direction) 
distribution to the true momentum (direction) distribu- 
tion using a Bayesian unfolding method [39( using the 
MC simulation to define the unfolding matrix. Figure [Hfl 
(|20p shows the true versus reconstructed ir° momentum 
(angle) distributions; these figures are graphical represen- 
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FIG. 16. Cosine of the reconstructed n angle with respect 
to the beam direction after all event selection cuts. 

tations of the smearing matrices used for the unfolding. 
Finally, we perform the efficiency correction to obtain the 
true 7T° momentum (direction) distribution. Figures [5Tj 
and [22] show the 7T° momentum and direction distribu- 
tions, respectively, after background subtractions, con- 
versions to the true ir° kinematics and efficiency correc- 
tions. To compare the shapes of the distributions, the to- 
tal numbers of entries in the distributions are normalized 
to unity both for the measurement and the MC expec- 
tation. The shapes of these two distributions agree with 
the MC expectation. The systematic errors of Fig.[2~Tland 
1221 are expected to arise from the same sources, and are 
estimated in the same manner, as described in Sec. IV A 4l 



FIG. 17. The true n momentum of generated MC events 
(solid) and selected MC events (dashed) for all NCw produc- 
tion processes (top), and the true 7r° momentum dependence 
of the efficiency of NCw production (bottom). 

C. Coherent pion production 

In the coherent pion production, the neutrino interacts 
the entire nucleus. In this case, the following relation 
should be satisfied. 

1*1 < ^ (8) 

where t and R are the momentum transfer to the nucleus 
and the nuclear radius, respectively. Using Eq. [51 we can 
deduce 

E„o{l - cos^o) < i ~ 100 MeV (9) 

according to Ref. [4(| ■ In this equation, the E„o and 9^o 
are the ir° energy and direction with respect to the neu- 
trino beam, respectively. Hence, the fraction of coherent 
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FIG. 18. The true -k angle with respect to the beam direc- 
tion for generated MC events (solid) and selected MC events 
(dashed) for NCw production (top) , and the efficiency for NC 
7T° production as a function of the true ty° direction (bottom). 
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FIG. 19. The true vr momentum (P true ) versus reconstructed 
7T° momentum (P rec ) from the MC simulation. The solid line 
shows the identity. 
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7T° production is extracted from the E™o(l — cosf5 rec ) 
distribution shown in Fig. [531 where E™a is the recon- 
structed 7T° energy calculated as E^f + E^f. We fit 
this distribution using three templates made by divid- 
ing the final MC sample into NC coherent 7T°, NC reso- 
nant 7r°and background samples. Two parameters, R co h, 
R ros , scale NC coherent 7r°and NC resonant 7r°templatcs 
independently. The scale of the background sample is 
fixed to unity. The expected number of events in the i- 
th (i = 1,2,.. .,N(= 20)) bin in the £^(1 - cosf5 roc ) 
distribution is expressed as: 



ivr p = r< 



3h x A^ coh + R rc 



x NJ CS + N, 



BG 



(10) 



FIG. 20. The true it direction (cos # true ) versus reconstructed 
7T° direction (cosf9 rec ) from the MC simulation. The solid line 
shows the identity. 



The fit minimizes the following x 2 '- 

f(N ohs ;N exp ) 
f(N ohs ; N ohs ) ' 



X 



-2 In 



(11) 



where 7V obs ( exp ) represents the observed 
(expected) number of events in all bins 

(iv ob s (cxp)^ob s (oxp) 5 _ _ 5Ar ob s (cxp) ) ^ ^obs.^exp) 

is the Poisson likelihood to find A^ obs events when 
jyoxp even ^ s are expected. When the systematic errors 
for each bin and their correlation expressed using the 
covariance matrix Vjk (j,k = 1,2, ...,N(= 20)) are 
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FIG. 21. The 7r° momentum distribution after all corrections 
described in the text, with statistical (error bars) and sys- 
tematic (red boxes) uncertainties. The dashed line shows the 
Monte Carlo expectation based on the Rein and Sehgal model. 
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variations drawn from Gaussian distributions about the 
expectations for each bin, with correlations, estimated 
from the MC simulation. Using Xi Jn for the m-th ex- 
pectation in the i-th bin, the likelihood is expressed as 



f(N ohs ;N^;V) = —J2l 



M ^ 11 iV obs 

m i 1 



(13) 



where M is the total number of random samples (1000). 
The result of the fit is: 



R coh = 0.56 ±0.34, 
R res = 1.33 ±0.16. 



(14) 
(15) 



The -E^5 C (1 — cos# rcc ) distribution after the fitting is 
shown in Fig. [23] The x 2 per degree of freedom (DOF), 
before the fit is 17.8/20 = 0.89, and it is 12.9/18 = 0.72 
after the fit. The statistical error and all systematic er- 
rors described in Sec. IV A 41 are included in the errors of 
R co h and R ros . Without the systematic errors, we ob- 
tain 0.79±0.30(stat.) and 1.24±0.13(stat.) for R coh and 
R rcs , respectively. The dominant systematic source is the 
detector response. 
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FIG. 22. The it angular distribution after all corrections de- 
scribed in the text, with statistical (error bars) and systematic 
(red boxes) uncertainties. The dashed line shows the Monte 
Carlo expectation based on the Rein and Sehgal model. 



given, the likelihood is expressed as 
f(N ohs ;N cxp ;V) = A 



~[dxi 

i=i 



_/V obs ! 



x exp 



1 



N N 



j=l k=l 



,(12) 



where A is the normalizing constant. To calculate this 
integral, we generate 1000 MC expectations with random 



FIG. 23. £; c c. c (l - cos(9 rc 
tribution and other NCw 
simulation. 



) after fitting. The coherent con- 
are separately shown for the MC 



The ratio of the NC coherent ir° production to the total 
CC cross sections from the MC prediction based on the 
Rein and Sehgal model is 1.21 x 10~ 2 . Hence, the cross 
section ratios are measured to be: 



a (NCcohTr ) 
ct(CC) 



= R coh x 1.21 x 10~ 2 , 
= (0.68 ±0.41) x 10~ 2 , 



(16) 



where R co h is 0.56±0.34. The mean neutrino energy for 
NC coherent ir° events in the sample is estimated to be 
1.0 GcV. This result is 1.6 standard deviations above the 
no coherent production assumption and consistent with 
the MC prediction based on the Rein and Sehgal model. 
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VI. CONCLUSIONS 

In conclusion, we have observed the production of the 
NC7T events by a muon neutrino beam on a polystyrene 
target (C 8 H 8 ) using the SciBooNE neutrino data set of 
0.99 x 10 20 protons on target. The ratio of the NCtt 
production to total CC cross sections is measured to 
be (7.7 ± 0.5(stat.) ± 0.5(sys.)) x 10~ 2 at mean neu- 
trino energy 1.1 GeV. The MC prediction based on the 
Rein and Sehgal model [i| is 6.8 x 10~ 2 . The measured 
shapes of the ir° momentum and angular distributions, 
as shown in Figures [21] and [22] agree with the MC pre- 
diction within uncertainties. The ratio of NC coherent 
7r°production to the total CC cross section is measured 
to be (0.7 ± 0.4 ) x 10~ 2 based on the Rein and Sehgal 
model [22!, while the MC prediction is 1.21 x 10~ 2 . 
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